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ABSTRACT 

Milli-second pulsars (MSPs) are rapidly spinning neutron stars, with spin periods 
P s < 10 ms, which have been most likely spun up after a phase of matter accre¬ 
tion from a companion star. In this work we present the results of the search for the 
companion stars of four binary milli-second pulsars, carried out with archival data 
from the Gemini South telescope. Based upon a very good positional coincidence 
with the pulsar radio coordinates, we likely identified the companion stars to three 
MSPs, namely PSR J0614-3329 (g=21.95±0.05), J1231-1411 (g=25.40±0.23), and 
J2017+0603 (g=24.72±0.28). For the last pulsar (PSR J0613—0200) the identification 
was hampered by the presence of a bright star (g=16±0.03) at ~2" from the pulsar 
radio coordinates and we could only set 3er upper limits of g = 25.0, r — 24.3, and 
i = 24.2 on the magnitudes of its companion star. The candidate companion stars to 
PSR J0614—3329, J1231—1411, and J2017+0603 can be tentatively identified as He 
white dwarfs (WDs) on the basis of their optical colours and brightness and the com¬ 
parison with stellar model tracks. From the comparison of our multi-band photometry 
with stellar model tracks we also obtained possible ranges on the mass, temperature, 
and gravity of the candidate WD companions to these three MSPs. Optical spec¬ 
troscopy observations are needed to confirm their possible classification as He WDs 
and accurately measure their stellar parameters. 
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1 INTRODUCTION 

Radio pulsars are interpreted as rapidly spinning and 
strongly magnetised neutron stars powered by their rota¬ 
tional energy. The most recent compilation (1.52) of the 
Australia National Telescope Facility (ATNF) Pulsar Cat¬ 
alogue (Manchester et al. 2005) lists over 2300 radio pul¬ 
sars. An important subgroup (~ 270 objects) is that of 
the so-called milli-second pulsars (MSPs), characterised by 
very short spin periods ( P B < 10 ms), high spin stability 
(P a ~ 10~ 18 10 -21 s s _1 ), spin-down ages P S /2P S ~ 1 10 
Gyrs, surface magnetic fields B ~ 10 8 -10 9 G, and spin-down 
energy E ~ 10 32 -10 36 erg s -1 . 

A large fraction of MSPs (164) are in binary systems 
and are mostly located in the Galactic plane. The fact that 
many MSPs are found in binary systems supports the com¬ 
monly accepted scenario for the formation and evolution of 
MSPs (the “canonical recycling scenario Alpar et al. 1982; 
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Bhattacharya & van den Heuvel 1991). This suggests that 
MSPs form in binary systems from old (slowly rotating) 
pulsars and that, at a later stage, they are reaccelerated 
to milli-second spin periods as the result of mass accretion 
from their non-degenerate companion stars. The pathways 
to the MSPs are different according to the initial mass ratio 
and orbital separation of the components, which drives the 
length of the accretion processes (see, e.g. Tauris 2011, 2015; 
Tauris et al. 2012). 

The most common pathway leads to a companion star 
that has almost or completely lost its external layers (white 
dwarf; WD) orbiting a rapidly spinning pulsar (Lyne et al. 
1987; Alpar et al. 1982; Bhattacharya & van den Heuvel 
1991; Lorimer & Kramer 2005). In this case, the companion 
star is expected to be either an He WD of mass 0.1 Mq < 
Me < 0.5A/q or a more massive Carbon-Oxygen (CO) WD, 
of mass 0.5Mq < Me < 1 Mq. When the companion star has 
not reached the WD stage yet it may be ablated by irradia¬ 
tion from the pulsar relativistic wind, which possibly leads 
to the formation of a solitary MSP. Two distinct families of 
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binary MSPs with non-degenerate companions emerge ac¬ 
cording to the degree of the ablation processes (e.g., Roberts 
2013). The first one is that of the so-called black widow 
(BW) MSPs, where the companion is a very low-mass star of 
Me < O.IMq, almost fully ablated by the pulsar wind. The 
second one is that of the redbacks (RB) MSPs, where the 
companion is only partially ablated and has an higher mass 
of Me ~ 0.1-0.4Mq. Therefore, the study of MSPs in binary 
systems is fundamental to understand the final stages of the 
binary pulsar evolution, including the different outcomes of 
the recycling process (e.g., Possenti & Burgay 2008; Possenti 
2013), and investigate the possible evolutionary connections 
between different MSP types (e.g., Benvenuto et al. 2014). 

In all these types of studies, the optical identification of 
the MSP companion star plays a crucial role (e.g., Pallanca 
et al. 2012, 2013, 2014; Mucciarelli et al. 2013). In particular, 
optical observations, either via spectroscopy or broad-band 
photometry, are key to determine the star classification and 
physical parameters, such as surface temperature, surface 
gravity, and chemical composition of the atmosphere. This 
information then allow one to determine the age of the bi¬ 
nary system (and that of the MSP independently of the mea¬ 
surements on P s and P s ) and track the system evolutionary 
history. Optical observations also allow one to find evidence 
for irradiation of the companion star through the identifica¬ 
tion of hot spots on the stellar surface. This can be a further 
tracer of the ablation process in BW and RB systems, af¬ 
ter the more direct one which comes from the observation 
of eclipses of the radio signal as it propagates through the 
stellar wind from the irradiated MSP companion. For only 
~ 50 binary MSPs the companion star has been identified 
in the optical. Thus, for only about one third of the MSPs 
we know the optical characteristics of the companion star. 
In most cases, this is mainly due to the relatively large dis¬ 
tances of these systems and the interstellar extinction in the 
Galactic plane, which hampers deep multi-band optical ob¬ 
servations. For the vast majority of non-eclipsing MSPs, the 
companion stars are found to be He WDs, although more 
cases of MSPs with CO WD companions have been found 
in the last few years (e.g., Mignani et al. 2014). 

In this manuscript we report on the search for the com¬ 
panion stars of four unidentified binary MSPs carried out 
using archival data from the Gemini-South Telescope. We 
describe the MSPs sample and the observations in Section 
2, while we present and discuss the results in Section 3 and 
4, respectively. 

2 OBSERVATIONS 
2.1 Target Description 

The names, measured radio coordinates, spin period (Ps), 
period derivative ( P s ) and inferred characteristics, spin- 
down age (r), dipolar magnetic field ( B ), spin-down energy 
( E ), of the four MSPs studied in this work are summarised 
in Table 1. 

Besides the radio band, all these MSPs have been de¬ 
tected as 7 -ray pulsars by Fermi (Abdo et al. 2013). PSR 
J0614—3329, J1231—1411, and J2017+0603 were discov¬ 
ered during radio follow-ups of unidentified Fermi sources 
(Ransom et al. 2011; Cognard et al. 2011), whereas PSR 


J0613—0200 was discovered in radio by Manchester et al. 
(1996), prior to its detection as a 7 -ray pulsar (Abdo et al. 
2009). 

All of them are also detected in the X rays by either 
Suzaku, Chandra, or XMM-Newton (see Abdo et al. 2013 
and references therein). All these pulsars, with the excep¬ 
tion of PSR J0614—3329, are in tight binary systems with 
orbital periods of ~ 2 days or less and have almost cir¬ 
cular orbits (Table 2), as expected from the outcome of 
the mass-accretion phase onto the MSP. The pulsar mass 
functions suggest low-mass companion stars with minimum 
masses Me ~ 0.13-0.3Mq, after assuming a pulsar mass 
Mp=1.4Mq and an orbital inclination i = 90° as indica¬ 
tive values. These assumptions are required since for none 
of these MSPs it was possible to measure the post-Keplerian 
parameters so far and for none of them evidence of radio 
eclipses has been found from the available observations (Ver- 
biest et al. 2009; Cognard et al. 2011; Ransom et al. 2011). 

Deep optical investigations with 8 m-class telescopes 
have never been reported for these four MSPs. For both 
PSR J0613—0200 and PSR J1231—1411, optical observa¬ 
tions were performed in January 2010 with the 2.4 m Isaac 
Newton Telescope (INT) at the La Palma Observatory (Ca¬ 
nary Islands) soon after their detection as 7 -ray pulsars 
(Collins et al. 2011), whereas for PSR J2017+0603, the 
only optical observations are those performed with the Swift 
Ultraviolet Optical Telescope (Cognard et al. 2011). In all 
cases, no candidate companion star to the MSP was de¬ 
tected. 

2.2 Observation Description 

We downloaded broad-band images of the MSP fields from 
the public Gemini science archive 1 . Observations were taken 
between August 5 2010 and January 29 2011 with the 
Gemini-South telescope on Cherro Pachon (Chile) under 
programme GS-2010B-Q-56. Observations were performed 
using the Gemini Multi-Object Spectrograph (CMOS). At 
the time of the observations, the camera was still mount¬ 
ing the original three-chip EEV CCD detector (2048x4068 
pixels each). This had a combined (unvignetted) field-of- 
view of 5l5 x 5f5, with gaps of 2 , /8 between each chip. 
The pixel scale was 0 / . / 1454 for a 2 x 2 binning. Observa¬ 
tions were performed through the g_G0325 (A = 4750A; 
AA = 1540A), r_G0326 (A = 6300A; AA = 1360A), and 
LG0327 (A = 7800A; AA = 1440A) filters, very similar to 
the g’, r’, and i’ used by the Sloan Digitised Sky Survey 
(SDSS; Fukugita et al. 1996). 

For all targets the observation sequences in each hlter 
consist of both short (typically 30 s) and long (240 s) expo¬ 
sures, with the former taken not to saturate bright stars in 
the field that can be used as secondary photometric refer¬ 
ence frames. The exposure sequence was performed applying 
a dithering of a few arc seconds along the X axis of the detec¬ 
tor both to compensate for the gap between the CCD chips 
and account for the fringing affecting the i band. All targets 
were observed close to the zenith, with airmass mostly be¬ 
low 1 . 2 , in dark time and mostly in photometric conditions. 
The complete observation log is reported in Table 3. 

1 http://cadcwww.dao.nrc.ca/gsa/ 
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Table 1. Coordinates, spin period P s and period derivative P s of the observed MSPs, together with the inferred values of the spin-down 
age (r), dipolar magnetic field ( B ) and rotational energy loss ( E ), as derived from the ATNF pulsar data base. For the MSP coordinates, 
the number in parentheses represent the la uncertainty on the last quoted digit. For those pulsars for which a proper motion has been 
measured (PSR J0613—0200 and J1231—1411), the values of P s , r, B and E , have been corrected for the Shlowskhi effect. 


Pulsar 

002000 
( hms ) 

^J2000 
(o ' ”) 

Ps 

(ms) 

Ps 

(10 18 s s -1 ) 

T 

(10 9 yr) 

B 

(10 8 G) 

E 

(10 34 erg cm 2 s 1 ) 

J0613—0200 1 

06 13 43.97514(1.1) 

-02 00 47.1737(4) 

3.06 

0.00882 

5.50 

1.66 

1.3 

J0614—3329 

06 14 10.3478(3) 

-33 29 54.118(3) 

3.14 

0.0175 

2.84 

2.38 

2.2 

J1231-1411 2 

12 31 11.3132(7) 

-14 11 43.63(2) 

3.68 

0.0228 

2.56 

2.93 

1.8 

J2017+0603 

20 17 22.7044(1) 

+06 03 05.569(4) 

2.89 

0.0083 

5.53 

1.57 

1.3 


1 fi a cos(5) = 1.84 =b 0.08 mas yr - 1 ; = —10.6 =b 0.2 mas yr 1 (Verbiest et al. 2009) 

2 fi a cos(8) = —100 ± 2 mas yr - fis = —30 d= 4 mas yr —1 (Ransom et al. 2011) 


Table 2. Orbital parameters (orbital period P^ and eccentricity e) of the pulsars listed in Table 1 as derived from the ATNF pulsar 
data base (second and third columns), together with the recomputed lower limits on the companion mass (Me) inferred from the system 
mass function and assumed a pulsar mass Mp=1.4M0 and an orbital inclination i = 90°. Proposed classification of the companion star 
are given in the following column. References for the binary system mass function measurements and proposed classifications are given 
in the last column. 


Pulsar 

Pb 

(d) 

e 

M C 

(M q ) 

Class. 

Refs. 

J0613—0200 

1.198 

0.0000055 

0.13 

He WD 

Verbiest et al. (2009) 

J0614—3329 

53.584 

0.0001801 

0.28 

WD 

Ransom et al. (2011) 

J1231—1411 

1.860 

0.000004 

0.19 

WD 

Ransom et al. (2011) 

J2017+0603 

2.198 

0.0000005 

0.18 

He WD 

Cognard et al. (2011) 


Table 3. Observation log for the Gemini MSP observations 


PSR 

Date 

yyyy-mm-dd 

Filter 

Number of Exposures 
x Exposure Time (s) 

Airmass 

seeing 

n 

J0613—0200 

2010-11-06 

g 

1x30 

1.14 

0.77 



r 

1x5+1x15+1x30 

1.14- 1.15 

0.62 



i 

1x5+1x30 

1.5 

0.61 


2010-12-29 

g 

7x30+7x60+7x120 

1.14- 1.18 

0.94- 1.20 



r 

7x30+7x60+7x120 

1.13-1.14 

0.81- 1.18 



i 

7x30+7x60+7x120 

1.14- 1.17 

0.63-1.00 

J0614—3329 

2010-09-08 

g 

10x240 

1.1-1.58 

0.72- 1.18 



r 

10x240 

1.09-1.50 

0.59-0.94 



i 

8x240 

1.12- 1.45 

0.69-0.83 


2010-09-14 

i 

1x30 

1.02 

0.85 



i 

2x240 

1.17- 1.19 

0.63 


2010-10-28 

g 

1x30 

1.02 

0.97 



r 

1x30 

1.02 

0.88 

J1231—1411 

2010-12-31 

g 

1x30 

1.25 

1.15 



r 

1x30 

1.24 

1.13 



i 

1x30 

1.24 

1.03 


2011-01-03 

g 

6x240 

1.33-1.37 

0.87- 1.05 



r 

1x96+4x240 

1.30-1.68 

0.80-1.18 



i 

4x240 

1.39-1.57 

0.76-0.81 


2011-01-29 

g 

6x240 

1.04- 1.08 

0.84- 1.21 



r 

6x240 

1.04- 1.1 

0.89-1.28 



i 

6x240 

1.04-1.1 

0.82- 1.13 

J2017+0603 

2010-08-05 

g 

1x92+1x200+8x240 

1.24- 1.35 

0.85-1.15 



r 

1x87+2x240 

1.24 

0.85 


2010-08-10 

r 

7x240 

1.49-1.66 

0.85-1.05 



i 

10x240 

1.33-1.47 

0.60-0.89 


2010-08-31 

g 

1x30 

1.24 

0.95 



r 

1x30 

1.24 

0.82 



i 

1x30 

1.24 

0.80 
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(c) J1231—1411 (d) J2017+0603 

Figure 1. CMOS r-band images of the four MSP fields. The circles correspond to the computed MSP radio positions. Since the error 
on the radio coordinates is below 1 mas (see Table 1) the uncertainty on the computed radio position is dominated by the accuracy our 
astrometry calibration (0”2; Sectn. 3.2). The circle radius corresponds to three times such an uncertainty for a better visualisation. For 
the other pulsars, the companion stars are marked by the ticks and labelled by the pulsar names. PSR J0614—3329 is located a few 
arcsec north of the disk of the bright galaxy ESO 365-1 (Lauberts & Valentijn 1989). All images have been obtained from the co-addition 
of all r-band science frames, apart from the PSR J0613—0200 one which has been obtained from the co-addition of short exposures only 
to avoid saturation of the bright star near the pulsar position. 


2.3 Data Reduction and Calibration 

We reduced the CMOS images using the dedicated GMOS 
image reduction package available in the Image Reduction 
and Analysis Facility (IRAF) package 2 . After downloading 
the closest-in-time bias and sky flat field frames from the 
Gemini science archive, we used the tasks gbias and gif lat 
to process and combine the bias and flat-field frames, respec¬ 
tively. We then reduced the single science frames using the 
task gireduce for bias subtraction, overscan correction, im¬ 
age trimming and flat field normalisation. These tasks work 
on a chip-by-chip basis and produce fully reduced images for 
each of the three GMOS CCDs. From the reduced science 


2 IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


images, we produced a mosaic of the three GMOS CCDs 
using the task gmosaic. Finally, we aligned the reduced im¬ 
age mosaics to compensate for the dithering pattern applied 
to the exposure sequence and average-stacked the aligned 
images with the task imcoadd to filter out cosmic ray hits. 
Unfortunately, the dithering pattern chosen for the observa¬ 
tions happened not to be adequate to account for the effects 
of the fringing in the i band. For this reason, the i-band im¬ 
ages obtained after frame stacking were still significantly af¬ 
fected by strong fringing. Applying the stacking to different 
image sub-groups did not improve the results significantly. 
To remove the fringing, we followed the recipe suggested by 
the GMOS science operation team 3 . We downloaded from 
the Gemini science archive the template i-band fringe images 


3 http: //ww. gemini. edu/sciops/instruments/gmos/calibration/ 
example-cal-data?q=node/10456 
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closest in time to the epochs of our observations. Then, af¬ 
ter flux normalisation, we subtracted the fringe images from 
the single science images using the IRAF task girmfringe. 
Finally, we applied the same procedures as above to mosaic 
and stack the de-fringed i-band science images. 

We computed the astrometry calibration of the optical 
images (both image mosaics and stacks) using the wcstools 4 
suite of programs that automatically match the sky coor¬ 
dinates of stars selected from, e.g the Two Micron All Sky 
Survey (2MASS; Skrutskie et al. 2006) or the Guide Star 
Catalogue 2.3 (GSC2; Lasker et al. 2008) catalogues with 
their pixel coordinates computed by Sextractor (Bertin & 
Arnouts 1996). After iterating the matching process and ap¬ 
plying a sigma-clipping selection to filter out obvious mis¬ 
matches, high-proper motion stars, and false detections, a 
pixel-to-sky coordinate transformation was computed us¬ 
ing a second order polynomial function and we obtained, for 
the ground-based images mean residuals of ~ O'! 2 in the ra¬ 
dial direction, using 50 bright, but non-saturated, 2MASS 
stars. To this value we added in quadrature the uncertainty 
atr = OH 08 of the image registr ation o n the 2MASS reference 
frame. This is given by at r =\Jn/Nsas (e.g., Lattanzi et al. 
1997), where Ns is the number of stars used to compute the 
astrometric solution, n =5 is the number of free parameters 
in the sky-to-image transformation model, as ~ O'! 2 is the 
mean absolute position error of 2MASS for stars in the mag¬ 
nitude range 15.5 ^ K ^ 13 (Skrutskie et al. 2006). After 
accounting for the O'/OIS uncertainty on the link of 2MASS 
to the International Celestial Reference Frame (Skrutskie et 
al. 2006), we ended up with an overall accuracy of ~0 / . , 2 
on the absolute optical astrometry, i. e. comparable to the 
CMOS pixel size. Given the exquisite accuracy of MSP ra¬ 
dio coordinates (Table 1) the uncertainty on the absolute 
astrometry calibration of the CMOS images obviously dom¬ 
inates the accuracy on the positional coincidence with a po¬ 
tential companion star. 

We computed the photometry calibration using as a 
reference sets of secondary photometric standard stars ex¬ 
tracted from the American Association of Variable Stars Ob¬ 
servers (AAVSO) Photometric All-Sky Survey 5 (APASS), 
directly identified in each CMOS image of our MSP fields. 
We preferred to follow this approach rather than computing 
the photometry calibration using images of standard star 
fields (Smith et al. 2007) because there were no standard 
star observations taken on the same nights as the science 
observations. Moreover, a direct on-the-frame calibration 
obviously allows one to compensate for variations in the 
sky transparency during the night. The APASS photometric 
survey (Henden & Munari 2014) is calibrated on the SDSS 
photometric system and is, then, well suited to calibrate our 
CMOS g, r, i-band images. Like the SDSS, our photometry 
is, then, in the AB system (Oke 1974). 


4 http://tdc-www.harvard,edu/wcstools/ 

5 http://www.aavso.org/apass 


3 RESULTS 
3.1 Astrometry 

By using the astrometric solutions obtained as described in 
the section above, we checked all co-added science frames 
in each filter to look for objects detected at the computed 
radio position of our target MSPs (Table 1), which can be 
be regarded as potential companion stars. For those pulsars 
which have a measured proper motion, i.e. PSR J0613—0200 
(Verbiest et al. 2009) and PSR J1231 —1411 (Ransom et 
al. 2011), we corrected the reference radio positions to the 
epoch of the Gemini observations. 

For three pulsars (PSR J0614—3329, -J1231 —1411, 
J2017+0603), we clearly detected an object at the radio co¬ 
ordinates (see Fig. 1; top right to bottom left). In all cases 
the match was very good, with the object centroid almost 
perfectly coincident with the computed pulsar radio posi¬ 
tion. In particular, the measured angular separation between 
the MSP coordinates and those of its candidate companion 
was only ~ 0”2-0”3, i.e. comparable with the astrometric 
uncertainty of the CMOS images (~ 0”2; Section 2.3). Thus, 
in each case we considered the association with the MSP 
safe. To formally quantify the goodness of the associations 
between the MSPs and their candidate companions, we com¬ 
puted the probability that they are due to a chance coinci¬ 
dence. We computed this probability as P — 1—exp(— npr 2 ), 
where r is the measured angular separation between the 
MSP coordinates and those of its candidate companion (in 
arc second), and p is the density of stellar objects around 
the MSP positions per square arc second. Since, in princi¬ 
ple, we do not know a priori the expected brightness range 
for the companion stars, we computed p without applying a 
selection in magnitude. For PSR J0614—3330 we computed 
a stellar density p ~ 0.028 and estimated a chance coinci¬ 
dence probability P ~ 0.0035, whereas for PSR J1231 —1411 
and PSR J2017+0603 we estimated P ~ 0.004 (p ~ 0.014) 
and P ~ 0.010 (p ~ 0.038), respectively. In all cases, the low 
chance coincidence probability would rule out spurious asso¬ 
ciations and would suggest that the candidate companions 
detected in the CMOS images are indeed associated with 
the MSPs. 

For the fourth pulsar (PSR J0613—0200), the radio po¬ 
sition fell slightly off-centre (~ 2" to the North) with respect 
to the centroid of the closest object detected in the CMOS 
images (Fig. 1; top left), although it is still within its PSF. 
However, the measured position offset is much larger than 
the 3a uncertainty on our computed astrometric solution, 
with the uncertainty on the pulsar radio coordinates being 
negligible (Table 1). Thus, based on our astrometric solu¬ 
tion, we infer that this object is not associated with PSR 
J0613—0200. Presumably, this object is the ’’brightish star 
nearby” the pulsar, mentioned in Table 1 of van Kerkwijk 
et al. (2005). This star was also identified by Mignani et 
al. (2014) in near-ultraviolet and optical images from the 
Galaxy Evolution Explorer ( GALEX ), the XMM-Newton 
Optical Monitor (OM), and the Swift UltraViolet and Opti¬ 
cal Telescope (UVOT) and ruled out as a candidate compan¬ 
ion to the pulsar on the basis of its positional offset from the 
radio coordinates. The same star was also detected in images 
taken with the INT (Collins et al. 2011) but was considered 
unassociated to PSR J0613—0200 for the very same reason 
as above. The relatively high chance coincidence probability 
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with the radio position also suggests that this star is unre¬ 
lated to the pulsar. By assuming a radius r ~ 2", i.e. equal 
to the angular separation between the star and the MSP ra¬ 
dio coordinates, we obtain P ~ 0.13 (p ~ 0.038). To directly 
verify a possible, though unlikely, association between this 
star and the pulsar, we compared their measured proper mo¬ 
tions. According to UCAC-4 (Zacharias et al. 2013), this star 
has a proper motion of ^i a cos(S ) = 7.2 ± 2.8 mas yr~ x and 
Pi = —4.6±2.9 mas yr -1 in right ascension and declination, 
respectively, which is marginally consistent with the pulsar 
radio proper motion, fj, a cos(5) = 1.84 ± 0.08 mas yr ^ 1 and 
Pi = —10.6 ± 0.2 mas yr -1 (Verbiest et al. 2009). However, 
the UCAC-4 measurement is still below the 3er significance 
so that the result of such a comparison is not conclusive. 
Thus, with all pieces of evidence pointing towards a chance 
coincidence, we conclude that the pulsar companion is most 
likely hidden in the PSF wings of this bright star. 


3.2 Photometry 

Since for all the target MSPs, the radio position is neatly 
located very close to the nominal CMOS pointing position 
at the centre of chip # 2 , we decided to use only images from 
the central chip for the following photometry analysis. This 
choice has the advantage of simplifying the procedure, with 
the tolerable cost of loosing the relatively small field of view 
of the two external chips, some of which are also affected 
by vignetting and partially occulted by the telescope auto- 
guider. Moreover, adopting a smaller, but still significantly 
large field of view, has the advantage of reducing the effects 
of possible differential extinction along the line of sight due 
to, e. g. clumps in the interstellar medium. 

We computed the photometry of the candidate MSP 
companion stars, as well as of the field stars, through a 
semi-automatic pipeline that makes use of the DAOPHOT 
II package (Stetson 1987, 1994) distributed as part of IRAF 
software environment. The pipeline consists of a set of 
scripts executing sequentially all the tasks needed to ob¬ 
tain a list of objects from a given image, with positions, 
magnitudes and errors. Namely, per each filter these tasks 
are: i) obtain a master list of detected objects through the 
analysis of the master image, chosen as the co-added and 
exposure-map corrected science frame to maximise the de¬ 
tection signal-to-noise, ii) obtain a reliable PSF model for 
all the images in the data set (single exposures and co-added 
science frames), iii) perform the PSF-htting photometry and 
aperture correction steps on both single exposures and co¬ 
added science frames to obtain magnitudes and errors of 
every object in the master list, iv) apply the photometric 
calibration as described in the previous Section, v) extract 
relevant data for the candidate MSP companion stars and 
some control objects selected among the closest neighbour 
stars, and, per each MSP field, vi) perform a variability anal¬ 
ysis of the different photometry measurements to spot pos¬ 
sible flux variations of the candidate MSP companion stars. 
The only step of our photometry pipeline which is run man¬ 
ually is the selection of the objects for building the PSF 
model in step ii), from the master list obtained from the 
master image. Once the starting sample has been selected, 
the actual PSF modelling is automatised. This procedure is 
a ’’custom” version of what more specialised programs do 


(see e.g. ALLFRAME; Stetson 1994) and has been repeat¬ 
edly tested in the past for similar works. 

The final result of our photometry pipeline is a cata¬ 
logue of positions (detector and celestial coordinates) and 
g, r, i flux measurements and errors for each MSP field 
and for each image (single exposures and co-added science 
frames) in our sample. Filter-to-filter matching was done 
in a trivial way because all the single-mosaic images are 
carefully registered onto the master. From this data set, we 
extracted time series for each candidate MSP companion 
(Section 3.3), identified from the astrometric matching with 
the radio positions, and built both colour-magnitude (CM) 
and colour-colour (CC) diagrams using mean magnitudes 
for both the candidate MSP companion stars and all stars 
detected around the MSP positions (Section 3.4). 

Using the same approach as above, we computed a 
model PSF for all the PSR J0613—0200 images (single mo¬ 
saics and co-added science frames). Then, we used it to sub¬ 
tract the star detected ~ 2 " from the pulsar radio position 
and search for possible objects hidden in its PSF wings. 
Unfortunately, the low signal-to-noise of the residuals after 
the subtraction of the star PSF did not allow us to find ev¬ 
idence of any excess of signal that could be associated with 
an object detected at the pulsar position. Thus, the actual 
companion to PSR .J0613—0200 would be undetected in the 
CMOS images. Only higher spatial resolution observations, 
either with Adaptive Optics or with the Hubble Space Tele¬ 
scope, would make its detection possible. After subtracting 
the star PSF we estimated the upper limits on the flux of the 
PSR J0613—0200 companion star using the co-added sci¬ 
ence images as a reference. The 3<r upper limits computed 
at the pulsar radio position are g = 25.0, r = 24.3, and 
i - 24.2. 

The companion stars to MSPs often feature optical vari¬ 
ability along the orbital period of the binary system owing, 
for instance to tidal distortion of the star (e.g. Orosz & van 
Kerkwijk 2003) or the formation of hot spots on the star sur¬ 
face irradiated by the pulsar (e.g., Stappers et al. 2001). For 
this reason, we searched for possible evidence of variability 
in the light curves of the candidate companion stars to PSR 
J0614—3329, J1231—1411, J2017+0603, obtained from the 
multi-epoch CMOS photometry (Table 3). As a criterion, 
we used only the photometric measurements obtained in 
those filters where the candidate companion star is brighter 
and is detected with the highest occurrence. For all MSP 
fields, these conditions are matched by the observations in 
the r filter. However, in all cases the reduced \r of the flux 
measurement distribution is quite low, with values of 0.18 
(J0614—3330), 0.24 (J1231-1411), 0.06 (J2017+0603), and 
is compatible with random variations only. Moreover, the 
small number of flux measurements, concentrated in short 
night fractions and repeated on two or three nights only, 
sometimes separated by up to a few weeks (see Table 3), 
largely under sample the known orbital periods. This is par¬ 
ticularly true in the case of PSR J0614—3329 which has 
the longest orbital period (Pb=53.584 d) among these three 
MSPs. Therefore, we cannot use the available multi-epoch 
flux information as a further piece of evidence to support the 
proposed identification between the MSPs and their candi¬ 
date companion stars, so far only based on the very good 
positional coincidence with the pulsar radio coordinates. 
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Table 4. Coordinates and multi-band photometry for the candidate companions to the three MSPs (column one), as measured in the 
GMOS images. Magnitude values were computed from the mean of all the available measurements and are in the AB system (Oke 1974). 
Column four gives the angular separation A r between the measured candidate companion and the MSP coordinates (Table 1). 


PSR 

a 

( hms ) 

8 

(O / «) 

A r 

g 

r 

i 

J0614—3329 

06 14 10.333 

-33 29 54.19 

0719 

21.95+0.05 

21.70+0.03 

21.58+0.03 

J1231—1411 

12 31 11.299 

-14 11 43.39 

0731 

25.40+0.23 

23.95+0.06 

23.35+0.11 

J2017+0603 

20 17 22.714 

+06 03 05.82 

0729 

24.72+0.28 

24.06+0.25 

23.84+0.17 


3.3 Colour and magnitude analysis 

Since no significant variations are seen in the objects’ 
fluxes, we computed mean optical magnitudes in the g, 
r, and i filters for the candidate companion stars to PSR 
J0614—3329, J1231—1411, J2017+0603, after applying a cr 
clipping algorithm to filter out measurements more strongly 
affected by night-to-night fluctuations. The mean mag¬ 
nitudes of the MSP candidate companion stars are (Ta¬ 
ble 4): g=21.95±0.05, r=21.70±0.03, i=21.58±0.03 (PSR 
J0614—3330), g=25.40±0.23, r=23.95±0.06, i=23.35±0.11 
(PSR J1231—1411), and g=24.72±0.28, r=24.06±0.25, 
i=23.84±0.17 (PSR J2017+0603). 

Incidentally, we note that the mean magnitudes of 
the star detected ~ 2" from the computed radio position 
of PSR J0613—0200 (Fig.l, top left) are: g=16.21±0.03, 
r=15.54±0.03, i=15.35±0.03. This suggests that, at the 
pulsar parallactic distance of 1.25±1.09 kpc (Verbiest et 
al. 2009), it would be an early main sequence (MS) star, 
whereas the minimum companion star mass and the or¬ 
bital parameters (Table 2) indicate a WD companion star 
to PSR J0613—0200, possibly an He WD (Lorimer et al. 
1995). Moreover, MSPs with possible early MS companions 
are extremely rare, with PSR J1903+0327 being the only 
certified case so far (Freire et al. 2011). However, this MSP 
has a quite large orbital eccentricity (e ~ 0.44) and a long 
orbital period [Pb= 95.17 d), whereas PSR J0613—0200 is in 
an almost circular orbit with a short period (Table 2). This 
indicates that the two binary systems followed different evo¬ 
lutionary paths, hence with different companions, with PSR 
J1903+0327 probably being part of a triple system in origin 
with both the MS star and a WD (Freire et al. 2011). There¬ 
fore, the colour and magnitude analysis give a further piece 
of evidence that this star is unrelated to PSR J0613—0200. 
An updated optical proper motion measurement for this star 
(see Sectn. 3.1), together with radial velocity measurements 
from optical spectroscopy and/or multi-epoch optical pho¬ 
tometry, would indisputably rule out its association with 
PSR J0613—0200. 

For the other pulsars (PSR J0614—3329, J1231—1411, 
J2017+0603) we used the mean magnitudes of the candi¬ 
date companion stars as a reference for their classification by 
analysing their locations in the observed (i.e. not corrected 
for the reddening) CM and CC diagrams. The observed CM 
and CC diagrams for the PSR J0614—3329, J1231—1411, 
and J2017+0603 fields are shown in Figure 2-4 (black filled 
circles). The locations of the candidate MSP companion 
stars are shown as blue filled circles with error bars. In order 
to reject outliers and include only high-confidence measure¬ 
ments, we plotted stars for which at least nine measurements 
per filter are available and with a < 0.08. 


Table 5. Pulsar dispersion measure (DM), obtained from the 
ATNF pulsar data base, and distance (D), inferred from the 
Galactic free electron density along the line of sight (Cordes & 
Lazio 2002), of the last three MSPs listed Table 1. For the dis¬ 
tance, we assumed an uncertainty of ± 20%. The fourth and fifth 
columns give the hydrogen column density Ah derived from the 
fits to the MSP X-ray spectra (Abdo et al. 2013) and the Galac¬ 
tic extinction in the pulsar direction (Ay) derived from the Ah 
using the relation of Predehl&Schmitt (1995). 


Pulsar 

DM 

(pc cm -3 ) 

D 

(kpc) 

A h 

(10 2 ° cm“ 2 ) 

Ay 

J0614-3329 1 

37.049 

1.88 

2.01 

u.oo_ 0 12 

11231-1411 

8.09 

0.43 

11.3 + 5.1 

0.63+0.28 

J2017+0603 

23.918 

1.57 

10 

0.56 


1 The distance based on the DM might be overestimated by a 
factor of two or more (Ransom et al. 2011). 


As seen, also owing to the limited star sample, all the 
observed CM diagrams show a significant scatter that makes 
it difficult to recognise the characteristic stellar sequences. 
This scatter is mainly due to the different distance of the 
field stars and, to a lesser extent, to their different redden¬ 
ing. Moreover, since the coordinates of these MSPs points 
at different directions in the Galaxy they sample different 
stellar populations, i.e. from the Galactic Disk, Bulge, and 
Halo. In the case of the PSR J2017+0603 field, however, 
the observed CM diagram shows evidence of a more defined 
stellar sequence. This suggests that most stars in the field of 
view belong to a more homogeneous stellar population, e.g. 
in an open cluster, although there is no known open cluster 
(or candidate) in the pulsar field. Alternatively, they might 
all be at a comparable distance from the Sun, as suggested 
by the line of sight to the pulsar field, which intercepts the 
Sagittarius spiral arm. For all pulsars, the location of the 
companion stars in the observed CC diagrams lies in most 
cases along, or close to the field stars, suggesting that the 
companions have no unusual colours. This might indicate 
that their colours would not be very much affected by an 
hypothetical irradiation from the MSP. 


4 DISCUSSION 

We compared the positions of the candidate companion stars 
in the observed CM and CC diagrams with stellar popula¬ 
tion models and template WD evolutionary tracks account¬ 
ing for both the distance to the pulsar and the interstellar 
extinction along the line of sight. For none of the three MSPs 
there is a distance measurement based on the radio paral¬ 
lax. Thus, with all due caveats, we used as a reference the 
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Figure 2. Observed CM (top row) and CC (bottom row) diagrams for all the stars in a ~ 2'0 x 5.'5 sky region around the PSR 
J0614—3329 position (filled black circles). The candidate companion star to the pulsar is marked by the black filled circle. Stellar 
sequences simulated from the Besangon models for different values of distance are shown in light and dark grey. In the CM diagrams 
the dark grey regions correspond to distance values within ±20% the assumed pulsar distance (see Table 5), whereas in the CC diagram 
they correspond to magnitudes within ± 0.05 the r-band magnitude of the pulsar candidate companion star. Theoretical evolutionary 
tracks for both He WDs (0.2, 0.25, 0.35 Mq; light blue dashed lines) and low-mass CO WDs (0.35, 0.4, 0AbM@; dark blue solid lines) are 
represented, computed from the models of Panei et al. (2007) for an age up to 4.8 Gyr and for the nominal values of the MSP distance. 
Masses increase from right to left. Evolutionary tracks for MS stars for both Z=0.02 (red solid lines) and Z=0.0001 (red dashed lines) 
are also shown for different mass ranges (0.5—1 Mq, with steps of 0.1 Mq) and for the nominal MSP distance. Masses decrease from 
top to bottom. Only in this case we assumed a distance value half of that obtained from the DM (see discussion in Ransom et al. 2011). 
Both the simulated stellar sequences from the Besangon models and the WD and MS evolutionary tracks were plotted assuming a null 
reddening. Reddening vectors are shown on the top of each panel. The lengths of the vectors correspond to the E(B-V) estimated from 
the hydrogen column density inferred from the fits to the MSP X-ray spectrum. We used the extinction coefficients of Fitzpatrick (1999). 


pulsar distance obtained from the dispersion measure (DM) 
and the Galactic free electron density n e along the line of 
sight (NE2001; Cordes & Lazio 2002). We assumed a re¬ 
alistic uncertainty of ±20% on the computed distance, as 
the authors recommended. This uncertainty is much larger 


than that derived from the errors on the DM but accounts, 
in most cases, for systematic effects related to the possi¬ 
ble under/over estimation of the Galactic electron density 
in certain directions. For PSR J0614— 3329, however, the 
uncertainty on the distance is probably much larger than 
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Figure 3. Same as Figure 2 but for the PSR J1231—1411 field. 


20 % since the line of sight to the pulsar is almost tangent to 
the Gum Nebula and the distance is likely half as estimated 
from the NE2001 model (Ransom et al. 2011). No direct 
measurement of the interstellar reddening along the line of 
sight is available for these pulsars. However, all of them have 
been detected in the X rays (Abdo et al. 2013). Therefore, as 
a reference, we estimated the reddening from the hydrogen 
column density IVh along the line of sight derived from the 
spectral fits to the X-ray spectra after applying the relation 
of Predehl & Schmitt (1995). Then, we computed the extinc¬ 
tion in the different filters using the extinction coefficients of 
Fitzpatrick (1999). We note that for both PSR J0614—3329 
and J2017+0603 the X-ray spectrum is poorly constrained 
(see Table 16 in Abdo et al. 2013), and so is the value of the 
column density Ah - In particular, for the PSR J2017+0603 
the Ah was set to the Galactic value in the pulsar direction 
and scaled for the pulsar distance (10 21 cm -2 ). We obtained 


a new estimate of the Ah from the pulsar DM by applying 
the linear fit between these two quantities computed by He 
et al. (2013). For the DM towards PSR J2017+0603 (23.918 
pc cm -3 ) the fit yields Ah = 7.17^2’ 3 5 x 10 2 ° cm -2 , where 
the errors are associated with the fit 90% confidence inter¬ 
val. The Ah value computed from the He et al. (2013) fit is 
consistent with that reported in Abdo et al. (2013) and cor¬ 
responds to Ay = 0.40^q'^ 2 ! after applying the relation of 
Predehl & Schmitt (1995). The DM, inferred distance, Ah, 
and estimated interstellar reddening along the line of sight 
for the three MSPs are summarised in Table 5. 

Firstly, we compared the observed CM and CC dia¬ 
grams of the MSP fields with simulated stellar sequences 
computed from the Besangon models (Robin et al. 2004). We 
simulated these sequences for different stellar populations, 
i.e belonging to the MS, Red Giant Branch (RGB), or the 
WD branch, and for distance values up to 15 kpc. The simu- 
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Figure 4. Same as Figure 2 but for the PSR J2017+0603 field. 


lated sequences for each MSP field are shown in the panels in 
Figure 2-4 as the grey scale maps. In the CM diagrams, the 
dark grey regions correspond to distance values within ±20% 
the assumed pulsar distance (see Table 5), whereas in the CC 
diagrams they correspond to magnitudes within ± 0.05 the 
r-band magnitude of the pulsar candidate companion star. 
As seen, the spread in the simulated stellar sequences well 
reproduces the spread in the observed points, as expected 
for different stellar populations at different distances. Since 
the field stars are affected by an unknown interstellar ex¬ 
tinction, and the simulations based on the Besangon models 
simply compute a reddening scaled proportionally to the as¬ 
sumed distance in a given direction, introducing a reddening 
correction in our simulations might bias a direct comparison 
between the observed and the simulated stellar sequences. 
Therefore, for simplicity, in all cases we simulated the stellar 
sequences assuming a null reddening. Each panel also shows 


the reddening vector corresponding to the E(B-V) estimated 
from the IVh measured along the line of sight to each MSP 
(Table 5) and computed using the extinction coefficients of 
Fitzpatrick (1999). Then, we used the reddening vectors as 
a reference to trace the extinction-corrected locations of the 
observed points for the MSP companion stars (blue points) 
along the simulated stellar sequences. 

As seen, for all the three MSPs the location of the can¬ 
didate counterpart in the diagrams falls off the region of the 
simulated MS and close to that of the simulated WD se¬ 
quence. This is consistent with a WD identification for all 
the candidate companion stars, as proposed in the literature 
(Ransom et al. 2011; Cognard et al. 2011). We note that non¬ 
degenerate companion stars as those in BW or RB systems 
would be much closer than observed to the late MS (see, e.g. 
Pallanca et al. 2012). Moreover, the minimum mass of the 
companion stars and the orbital period of the binary system 


MNRAS 000, 1-17 (2015) 














Gemini optical observations of binary millisecond-pulsars 11 


(Table 2) would rule out that these MSPs are BWs or RBs 
(see, e.g. Figure 1 in Roberts 2013). Indeed, the minimum 
companion masses are larger than those of BW companions, 
whereas the orbital periods are larger than those of known 
RB binary systems, which have all orbital periods shorter 
than 1 d. The latter case is especially true for PSR J0614- 
3329 which has an orbital period of 53.584 d. 

In order to better determine the WD characteristics, 
obtain a first tentative estimate of their mass and an inde¬ 
pendent estimate of the age of the binary system, we used 
the model evolutionary tracks of Panei et al. (2007), which 
are computed for both He and low-mass CO WD types, and 
different WD masses. In particular, these tracks are quite 
suited to our goal since the assumed range of WD masses is 
close to the minimum mass of the companions inferred from 
the mass function of the binary systems (Table 2). These 
tracks 6 also predict magnitude values for both WD types 
in the same photometric system as used by the G-MOS ob¬ 
servations, i.e. the SDSS one (Fukugita et al. 1996). As an 
example, Figure 2-4 show the evolutionary tracks simulated 
for both an He WD (light blue) and a low-mass CO WD 
(dark blue). 

The evolutionary tracks span an age range up to 4.8 
Gyr, where the age limit corresponds to that of the avail¬ 
able model (Panei et al. 2007), and different masses. As a 
reference, we also plotted the evolutionary tracks of MS stars 
for both Z=0.02 (red solid lines) and Z=0.0001 (red dashed 
lines) and for different mass ranges (0.5-1 Mg, with steps of 
0.1 Mg). In each figure, the lines correspond to the nominal 
value of the MSP distance (Table 5). Following Ransom et 
al. (2011), for PSR J0614—3329 we assumed a distance value 
half of that obtained from the DM and the NE2001 model. 
As done for the simulated stellar sequences computed from 
the Besangon models, we plotted the simulated WD and MS 
evolutionary tracks assuming a null extinction, for simplic¬ 
ity. As seen, for both PSR J0614—3329 and J2017+0603 the 
location of the candidate companion star in the CM and CC 
diagrams seems to be more compatible with the simulated 
evolutionary tracks for an He WD (light blue lines), than for 
a low-mass CO WD (dark blue lines). The difference between 
the He and CO WD evolutionary tracks is more evident for 
PSR J0614—3329, whereas for J2017+0603 the r-i colour is 
partially compatible with both WD types. The location of 
the candidate companion star to PSR J1231 —1411 in the 
CM and CC diagrams does not overlap with any of the evo¬ 
lutionary tracks, although it seems to be consistent with the 
extrapolation of the He WD evolutionary tracks. Since these 
are computed for a maximum age of 4.8 Gyr, it means that 
the companion star must be significantly older than those 
of PSR J0614—3329 and J2017+0603. To summarise, our 
qualitative analysis of the CM and CC diagrams seems to 
suggest that the candidate companions to these three MSPs 
are more probably He WD than low-mass CO WDs. How¬ 
ever, owing to the uncertainties in extrapolating the evolu¬ 
tionary tracks for ages above 4.8 Gyr, we regard the possible 
identification of the PSR J1231 —1411 companion as an He 
WD as more uncertain than in the other cases. In addition, 
it should be considered that CO WDs tracks tend to sep- 


6 http://fcaglp.fcaglp.unlp.edu.ar/evolgroup/TRACKS/ 
tracks Jieliumcore. html 


arate and spread out at low effective temperatures due to 
different atmopheric compositions. An example of this can 
be seen in Bergeron et al. (2011) and refs, therein. 

We tried to verify more quantitatively a possible as¬ 
sociation with either of the two types for different values 
of the WD mass and age, and for different values of the 
distance and reddening. To this aim, we systematically in¬ 
vestigated all possible combinations in a four dimensional 
parameter grid and selected only those for which the optical 
magnitudes and colours predicted by the model evolutionary 
tracks for He or low-mass CO WDs were found to be consis¬ 
tent with the observed ones, within the measured photom¬ 
etry uncertainties. We considered the case of He WDs first. 
Fig. 5, 6 and 7 show the allowed combinations for each of the 
selected pairs of parameters for the He WD case. For visu¬ 
alisation purposes, in each panel the regions corresponding 
to the allowed parameter combinations are plotted in light 
grey. The discontinuity in the grey regions noticed in some 
panels are purely an effect of the parameter quantisation. 
The black thick lines in each panel mark the investigated 
parameter range. We considered a range of values around 
those expected from radio and optical observations, with 
a generous tolerance to account for the associated uncer¬ 
tainties. In particular, we considered a range of distances 
0.01 < d < 2 kpc (0.01 kpc steps) and 0 < Ay < 1 (0.01 
magnitude steps). For the companion mass we considered 
the range O.188M 0 < Me < O.488M 0 (0.002 M 0 steps), 
which is that of the model evolutionary track for He WDs 
(Panei et al. 2007). We derived WD tracks for intermediate 
mass values from a linear interpolation between the tracks 
corresponding to the available mass values. Finally, for the 
age we considered a range of values obtained from the com¬ 
puted WD tracks, for a given value of WD mass and extinc¬ 
tion, with steps of 0.1 magnitudes in the r-band magnitude. 
We remind that for PSR J1231 —1411 the comparison with 
the WD tracks is based on their extrapolation for ages above 
4.8 Gyr and is obviously more uncertain than for the other 
two MSPs. 

Possible parameter ranges, purely imposed by the ob¬ 
served colours and magnitudes of the companion star, are 
represented by the blue dashed lines intercepting tangen¬ 
tially the light grey regions in each panel in Fig. 5-7. Note 
that in some panels the black thick lines and the dashed blue 
lines coincide. This corresponds to those cases where the 
uncertainty on our photometry is not sufficiently small to 
better constrain the parameter value. Then, we considered, 
with the due caveats, the limits imposed by the radio obser¬ 
vations on the MSP parameters, such as the lower limit on 
the companion mass Me (Table 2) and the pulsar distance 
(Table 5), and the limits imposed by the X-ray observations, 
such as the Ah and the inferred interstellar extinction (Ta¬ 
ble 5). We did not apply any limit to the WD age based on 
the pulsar spin-down age because of the well known difficul¬ 
ties in determining a reliable uncertainty range for this value 
other than that the formal error derived from the measured 
P s and P a , of course (e.g., Lorimer & Kramer 2005). Finally, 
we selected only the configurations for which all conditions 
imposed to all pair of parameters by the limits derived from 
optical, radio, and X-ray observations were simultaneously 
satisfied. The selected parameter configurations are repre¬ 
sented by the dark grey regions in the panels in Fig. 5-7. For 
PSR J1231-1411, the obtained age range originally spanned 
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Figure 5. Pair combinations of the four parameters which we used to fit the observed magnitudes and colours of the candidate companion 
star to PSR J0614-3329 with mode evolutionary tracks (Panei et al. 2007). This plot refers to the case of an He WD 

. The thick black lines mark the investigated region in each plane of the parameter space. The light grey regions mark, for each pair of 
parameters, the allowed parameter configurations which reproduce the observed optical properties of the candidate companion stars 
(magnitudes, colours). The extreme of these regions projected on each axis are marked by the blue lines, while the green dashed lines 
mark the parameter range constrained by the radio and X-ray measurements. Although the radio/X-ray constraints of each parameter 
pair seem to be more stringent than the optical ones, imposing that all conditions are simultaneously satisfied shrinks the allowed 
parameter configurations to the dark the grey regions. The dotted green line corresponds to the pulsar spin-down age. 


values larger than the age of the Universe. Therefore, we im¬ 
posed that the WD age is smaller than an arbitrary value, 
which we set to 10 Gyr. The plots in Fig. 6 were updated 
accordingly. 

We repeated the same analysis as above for the low- 
mass CO WD case. We assumed the same range of distance 
and extinction values as in the previous case and the range 
of companion star masses (O.351M0 < Me < 0.448 Mq) 
from the model evolutionary tracks for low-mass CO WDs 
(Panei et al. 2007). We computed the age range exactly as 


in the previous case. We found that all the allowed config¬ 
urations (not shown here) imply companion masses of <0.4 
M q (PSR J0614-3329, PSR .J2017 + 0603) and <0.36 M 0 
(PSR J1231-1411), i.e. at the low end of the mass range even 
for low-mass CO WDs. This would make a CO WD com¬ 
panion for the three MSPs somewhat less likely, although 
this possibility can be firmly ruled out only with a better 
characterisation of the star spectra through follow-up spec¬ 
troscopy observations. Our results are in line with the qual¬ 
itative analysis of the CM and CC diagrams (Figure 2-4) 
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Figure 6. Same as Fig. 5 but for PSR J1231-1411. 


and the possibility that the candidate companion stars of 
the three MSPs are He WDs. 

The derived parameter ranges for the three MSP com¬ 
panions are summarised in Table 6, for the He WD case. 
In the first half, we report the ranges of WD age, mass, dis¬ 
tance, and extinction obtained from the above analysis (Fig. 
5-7), whereas in the second half we report the values of the 
WD surface temperature, luminosity and surface gravity ex¬ 
tracted from the model, which correspond to the obtained 
age and mass ranges. Unfortunately, for both PSR J1231- 
1411 and PSR J2017+0603, which are the MSPs with the 
faintest companion stars (Table 4), the photometry errors 
were so large that several sets of parameters were able to 
fit the observed optical magnitudes and colours, although 
the parameter degeneracy has decreased (see Fig. 6 and 
7). On the other hand, the smaller photometry errors for 


PSR J0614-3329 made it easier to find possible ranges for 
the different parameters. We note that our results are influ¬ 
enced by the model uncertainties, including the characteri¬ 
sation of the WD atmosphere and composition, and by the 
overall uncertainties on the predicted magnitudes, which are 
more difficult to quantify than the photometry errors on the 
observed magnitudes, by some arbitrary assumptions on the 
parameter uncertainties, such as on the pulsar distance, and 
by the uncertainty on the extrapolated parameters, such as 
the interstellar extinction. Therefore, the inferred ranges for 
the WD parameters should not be taken rigidly but only as 
indicative. This is particularly true in the case of PSR J1231- 
1411 (see above). Spectroscopy observations and multi-band 
light curves will be needed to obtain a better characterisa¬ 
tion of the companion star properties and pave the way to 
a more robust comparison with stellar models. 
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Figure 7. Same as Fig. 5 but for PSR J2017 + 0603. 


For all our pulsars, we found a possible range of val¬ 
ues for the companion mass Me. As it has already been 
observed, a correlation has been found between the orbital 
period Pb and the companion mass Me for binary systems 
hosting a MSP and an He WD (see, e.g. Corongiu et al. 2012 
and references therein). Such a correlation was theoretically 
found by Tauris & Savonije (1999) and has been confirmed 
by the mass measurements of the companion star in about a 
dozen of such systems. The predictions by Tauris & Savonije 
(1999) indicate, for the three MSP binary systems discussed 
here, a value for the companion mass that would lie in the 
range that we obtained from our estimates (Table 6). This is 
seen in Fig. 8 (updated from Corongiu et al. 2012), where we 
show the Pb-Mc plot for all MSPs with an He WD compan¬ 
ion of measured mass together with the estimated masses for 


the companion stars to PSR J0614-3329, PSR J1231-1411, 
and PSR J2017 + 0603. 

For both PSR J0614-3329 and PSR J2017+0603, the in¬ 
ferred age estimate for the WD companion would be smaller 
than the pulsar spin-down age r by a factor of ~2.4-5 (Table 
1). On the other hand, in the case of PSR J1231-1411 the 
minimum age of the WD companion obtained from our anal¬ 
ysis would be ~ 70% greater than the pulsar spin-down age. 
However, it has been shown (e.g., Tauris 2012) that spin- 
down ages are not reliable age estimators for MSPs. There¬ 
fore, we must take this discrepancy with the due care. Since 
for PSR J1231-1411 we had to impose an arbitrary upper 
limit to the age of the companion star, only for PSR J0614- 
3329 and PSR J2017 + 0603 we could obtain a first estimate 
of the cooling age of the WD, which would be of the order 
of about 0.8 and 1.7 Gyr, respectively. 
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Orbital period vs companion mass for MSP—HeWD binaries 



M c (M 0 ) 

Figure 8. Orbital period F\, vs. companion mass Mq for MSP-Hc-WD binaries (green filled circles; updated from Corongiu et al. 
2012). Error bars correspond to 1 rj uncertainties in the masses. If not visible, their size is smaller than the symbol. The red diamonds 
correspond to the three MSPs for which we identified the companion stars: PSR J0614-3329, PSR J2017 + 0603, and PSR J1231-1411 
(top to bottom). The black filled triangle corresponds to the BW PSR B1957+20. The solid, dashed, and dot-dashed lines are the tracks 
corresponding to the theoretical correlation found by Tauris & Savonije (1999). 



PSR J0614—3329 

PSR J1231—1411 

PSR J2017+0603 

M C [M 0 ] 

0.276-0.414 

0.188-0.384 

0.188-0.42 

Age [Gyr] 

0.55-1.18 

4.30-10 

1.31-2.44 

D [kpc] 

0.76-1.02 

0.35-0.51 

1.25-1.87 

Ay [mag] 

0.24-0.71 

0.35-0.91 

0.2—0.4 

T [K] 

7237-9666 

2652-4009 

5412-6437 

log [ L/Lq] 

[-2.898, -2.386] 

[-4.549, -3.920] 

[-3.346, -3.014] 

log g 

7.107-7.572 

6.973-7.603 

6.778-7.625 


Table 6. Parameters for the candidate companion stars to the MSPs obtained from Fig. 5—7 and the comparison with He WD model 
evolutionary tracks (Panei et al. 2007). Top section: derived ranges for the mass and age of the companion star (Me) to the MSP, distance 
(D) and interstellar extinction along the line of sight (Ay). These ranges basically correspond to the outermost values of the dark grey 
regions of Figures 5, 6 and 7. In case of PSR J1231-1411, these ranges were obtained by imposing that the WD age is < 10 Gyr. Bottom 
section: ranges of WD surface temperature T, luminosity L, and gravity log g corresponding to the above mass and age ranges. 


For PSR J0614-3329 the derived range for the pulsar 
distance purely imposed by the observed colour and magni¬ 
tude of the companion star (blue lines in Fig. 5), would con¬ 
firm that the distance to the pulsar is indeed about half the 
value inferred from the DM and the NE2011 model (Table 
5), as suggested by Ransom et al. (2011). On the contrary, 
the corresponding range for PSR J1231-1411 would indicate 
that the pulsar distance cannot be much larger than the 
value obtained from the DM, whereas its high Galactic lati¬ 
tude (b ~ 48°) and the low value of the DM (8.09 pc cm -3 ) 
might have suggested that such a value is underestimated 
by a factor of two (Gaensler et al. 2008; Chatterjee et al. 
2009). 

We used the mass and age ranges of the WD compan¬ 
ions derived above to select the associated ranges of the star 
surface temperature T, luminosity and surface gravity log g 
(second half of Table 6) from the He WD model evolution¬ 
ary tracks (Panei et al. 2007). As seen, the WD companion 
star to PSR J1231-1411 would stand out for its lower tem¬ 
perature (T ~ 2600-4000 K), definitely on the lower end 
of the WD temperature distribution, and lower luminosity 
(log [L/Lq] « —4) with respect to the other MSP compan¬ 


ions. The inferred lower temperature was also anticipated 
by the observed reddish colour of the star (Fig. 3). Thus, 
the WD companion to PSR J1231-1411 would be one of the 
coolest WDs known so far. If confirmed, this result would 
be exceptional but not unheard of. We note that a very low 
temperature WD companion (T < 3000 K) has been likely 
identified also in the MSP binary system PSR J2222—0137 
(Kaplan et al. 2014). On the other hand, both the WD 
companion stars to PSR J0614-3329 and PSR J2017 + 0603 
would be hotter, with temperatures more compatible with 
the values expected for a WD (T ~ 7000-10000 K), and are 
more luminous ( log [L/Lq] « —3). 


5 SUMMARY AND CONCLUSIONS 

Using the Gemini-South telescope, we identified the 
likely companion stars to three MSPs (PSR J0614—3329, 
J1231 —1411, and J2017+0603). For a fourth pulsar (PSR 
J0613—0200), the identification of the companion star was 
hampered by the presence of a bright star (g=16±0.03) at 
~ 2" from the pulsar radio position. We ruled out with a 
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reasonable confidence that this star is a candidate compan¬ 
ion to the pulsar. The companion stars of the other three 
MSPs can be possibly identified as He WDs on the basis 
of multi-band photometry and comparison with simulated 
stellar models (Panei et al. 2007). For PSR J1231 —1411 the 
identification is somewhat more uncertain owing to the more 
ambiguous match with the evolutionary tracks. If these iden¬ 
tifications were confirmed, they would prove that He WDs 
tend to be the most common companions to non-eclipsing 
MSPs, as suggested by the number of MSP companions 
which have been identified in the optical so far. For PSR 
J0614—3329, J1231 —1411, and J2017+0603, we derive pos¬ 
sible ranges for the WD mass, age, surface temperature, lu¬ 
minosity, and gravity (see Table 6) within the ranges pre¬ 
dicted by the assumed He WD models and optical/radio 
observations. For none of these three MSPs we could look 
for flux modulations at the orbital period of the binary sys¬ 
tem, owing to the sparse data points and the uneven time 
coverage of the orbital period provided by the Gemini ob¬ 
servations. Future variability studies of the companion star 
fluxes through optical photometry, together with the mea¬ 
surements of the radial velocity curves through optical spec¬ 
troscopy with 8m-class telescopes, will firmly prove the pro¬ 
posed associations with the MSPs, currently based upon a 
very good positional coincidence with the radio coordinates 
and the colours of the candidate companions stars. Optical 
spectroscopy would also be crucial to verify our tentative 
classification of the MSP companion stars as He WDs, so 
far only based upon optical broad-band photometry in three 
bands, obtain more robust measurements of their stellar pa¬ 
rameters, such as the mass, surface temperature and gravity, 
and provide information on the atmosphere composition. 
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